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Abstract 



I discuss a prospect for mode identification from two-passband photometry 
of forthcoming BRITE space mission. Examples of photometric diagnostic dia- 
grams are shown for three types of main sequence pulsating variables: f3 Cephei, 
Slowly Pulsating B-type and 5 Scuti stars. I consider also taking into account 
. !^ the radial velocity data from simultaneous spectroscopy, which can be carried 
■ out from the ground. With such observations, much better discrimination of 
H I the spherical harmonic degree, £, can be accomplished and more constraints on 
- - ' stellar parameters and input physics can be derived. 



Introduction 

Nowadays space based observations allow detecting oscillation modes with lower 
and lower photometric amplitudes. We are already at the detecton threshold 
of the order of 10^^ — 10^^ rnag, resulting in a growing number of frequency 
peaks. There are many examples of excellent work both observational and 
theoretical based on WIRE and MOST data (e.g. Bruntt et al. 2007, Walker et 
al. 2005, Barban et al, 2007, Saio et al. 2007). Now we expect similar results 
from the just initiated COROT mission. 

However, for using these rich frequency data for asteroseismic modelling, 
mode identification is a prerequisite. In the case of main sequence pulsators, 
we are still far from obtaining very regular patterns in the oscillation spectra, 
which could help in solving the problem. We are also far from nonlinear theory, 
which would answer a question about mode selection mechanism. The main 
unresolved problem is why most of the theoretically unstable modes are not 
observed. Nowakowski (2005) suggested that the dominant effect of limiting 
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the mode amplitude is a collective saturation of the opacity driving mechanism, 
instead of a resonant mode coupling. 

In order to identify modes in (3 Cephei, 6 Scuti or Slowly Pulsating B- 
type stars, we need additional observables from multicolour photometry or/and 
spectroscopy. The photometric method of mode identification consists in using 
amplitude ratios and phase differences in different passbands; it is based on the 
semi-analytical formula for the light variation due to linear pulsation derived by 
Dziembowski (1977). Balona &l Stobie (1979) showed that modes with differ- 
ent values of ^ are located in different regions on the amplitude ratios vs phase 
differences diagram. Since then, the method has been applied to various types 
of pulsating variables by many authors (Watson 1988, Garrido et al. 1990, 
Heynderickx et al. 1994). The next important improvement was including 
nonadiabatic calculations by Cugier, Dziembowski & Pamyatnykh (1994), who 
applied the method to /3 Cep stars. Then, Balona &i Evers (1999) emphasized 
the problem of very high sensitivity of photometric amplitudes and phases to 
the treatment of convection in the case of 5 Set stars. A photometric identi- 
fication of t for SPB stars, based on nonadiabatic calculations, was performed 
by Townsend (2002). All these works were done by assuming that the rotation 
does not influence pulsation. However, main sequence pulsators are very often 
rapid rotators. The next step in developing the method was its extension to 
close frequency modes coupled by a fast rotation by Daszyriska-Daszkiewicz et 
al. (2002). The photometric method in the version for long-period g-modes 
in rotating stars, for which perturbation approach is no longer adequate, was 
formulated by Townsend (2003a) and Daszyriska-Daszkiewicz et al. (2007). 

A few years ago, Daszyriska-Daszkiewicz et al. (2003,2005a) proposed a 
new method of the identification of i, which uses the amplitudes and phases 
themselves and combines photometry, and radial velocity data. The method 
allows also to extract simultaneously a new asteroseismic probe, which yields 
constraints on stellar parameters and input physics, e.g., convection, opacities. 

BRITE (BRIght Target Explorer) is the first space-based mission which will 
perform two-colour photometry of bright stars. It will give an opportunity of not 
only detecting low-level oscillations but also of identifying their degrees I. How- 
ever, much more could be achieved, if ground-based spectroscopic observations 
are organized simultaneously. Then another type of information, contained in 
the radial velocity and line profile variations, would be supplied. 

The aim of this paper is to show what can be done for mode identifica- 
tion from observations with the two BRITE passbands. I will recall the basic 
formulae and show examples of photometric diagnostic diagrams for models of 
(3 Cephei, SPB and 5 Scuti stars. Then I will present diagrams which include 
amplitudes and phases from photometric and radial velocity variations. Finally, 
I will discuss uncertainties arising from effects of rotation, convection and at- 
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mospheric models. A summary and conclusions are given in the last section. 



Mode identification from photometry 

In order to calculate photometric amplitudes and phases, two inputs are needed. 
They come from 

• nonadiabtic theory of stellar pulsation, 

• models of stellar atmospheres. 

We assume linear pulsation theory, which is adequate because of small mode 
amplitudes in main sequence pulsators, and we use temporally static, plane par- 
allel atmosphere, which is justified because the eigenfunctions of the considered 
modes are nearly constant in the atmosphere. 

Let us consider a pulsation mode in the zero-rotation approximation, the 
geometry of which can be described by a single spherical harmonic, Y™, with 
the degree £ and the azimuthal order m. The shape of the radial eigenfunc- 
tions of the mode are determined by its radial order n. Then the local radial 
displacement is given by 

-=eYr{e,cl>)e-'^', (1) 

where e is the intrinsic mode amplitude, co is the angular pulsation frequency, 
which of course depends on {n,i,m); other symbols have their usual meanings. 
The corresponding changes of the bolometric flux, J-'ho\< and the local gravity, 
g, are given by 

^=£/rr(0,0)e-'-*, (2) 

>bol 

^ = -.(2 + ^)rr(M)o-. (3) 

The complex parameter, /, describes the ratio of the local flux perturbation to 
the radial displacement at the level of the photosphere and it is obtained from 
nonadiabatic calculations. 

The complex photometric amplitudes of the light variation in the passband 
A due to a pulsation mode with frequency oj can be written in the following 
form 

Axii) = -imGeVrii, Q)b}{Dl, + D^^e + Dl,) (4) 

where 

, _ 1 d\og{TAb^,\) 
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D2,e = {2 + £){!-£), (56) 
01og(-FA|b,"|) ... 

and z is the inclination angle. The partial derivatives of log{ J^x\b^\) over effec- 
tive temperature and gravity are derived from atmospheric models and 



(6) 



is the disc-averaging factor, containing the information about the visibility of 
the mode with a given degree i. The integrals are weighted by the limb- 
darkening law, h\{ii). The term D^^ describes the temperature effects, the 
term D2.1 stands for the geometrical effects, and the influence of gravity changes 
is contained in the term D^^. The terms ^ and ^ include the perturbation 
of the limb-darkening, and their ^-dependence arises from the nonlinearity of 
the limb-darkening law. 
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Figure 1: The disc-averaging factor, b(, as a function of the spherical harmonic degree, 
for the BR Johnson filters. The adopted stellar parameters are logTeff = 3.875 
and logs' = 4.0. 

The amplitudes and phases of the light variation are given by and 
a,Yg{A\), respectively. Having these numbers, we can construct photometric 
diagnostic diagrams in the form A^/Ay vs. (px — fy, where x and y denote 
passbands. These observables are independent of the intrinsic amplitude, e, 
inclination angle, i, and azimuthal order, m, because the product £^7" drops 
out. These is an advantage for the identification of t but also a disadvantage 
because the order m is beyond of the reach of the photometric method. 
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In this section I will give examples of photometric diagrams in the BR 
Johnson filters which are not very different from the BRITE passbands: BTl 
(390-460 nm) and BT2 (550-700 nm). The well known property of the factor 
6^ is that it decreases very rapidly with growing degree £; this can be seen 
from Fig. 1. In this paper I will consider modes with degrees up to £ = 6. 
All calculations were done with the Warsaw-New Jersey evolutionary code and 
nonadiabatic pulsation code of Dziembowski (1977). I used OPAL opacity ta- 
bles of Iglesias &L Rogers (1996) and the solar chemical composition of Grevesse 
&L Noels (1993), assuming the metallicity Z — 0.02. I adopted Kurucz atmo- 
spheric models in the NOVER-ODFNEW version (Castelli &i Kurucz 2004), 
which have more smooth flux derivatives than the standard Kurucz models, 
and the Claret nonlinear formula for the limb-darkening law. The standard 
atmospheric metallicity, [m/H]=0.0, and the microturbulence velocity, = 2 
km/s, were assumed. 

/? Cephei pulsators 
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Figure 2: The locations of unstable modes with degrees up to 6 for (3 Cephei star 
models of 12 A-Iq on the diagnostic diagrams involving Johnson B and R filters. The 
left panel contains all models in the main sequence evolutionary phase, and the right 
panel, the model with logTcff = 4.400 and logg = 3.89. 

I considered stellar models with a mass of M = 12Mq during main sequence 
phase of evolution, corresponding to the temperature range of log Toh = 4.445— 
4.347. In Fig. 2, the photometric diagram in the BR passband is presented. 
The left panel shows all modes which become unstable between ZAMS and 
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TAMS and the right panel shows unstable modes for only one stellar model 
with logTcff = 4.400 and \ogg — 3.89. Modes with different degrees £ are 
located in separated regions. The radial modes are spread over a wide range 
of the amplitude ratios and phase differences, whereas the nonradial modes are 
concentrated in small areas. This behaviour results from different contributions 
of the temperature (D^^/) and the geometrical effects {D2/) to the light 
variation (Daszyriska-Daszkiewicz et al. 2002) 

Slowly Pulsating B-type stars 

In order to calculate SPB oscillation, I chose the models with M = BMq on 
the main sequence, which include log Toff from 4.235 to 4.134. Positions of 
modes with different values of ^ on the photometric BR diagram are shown in 
Fig. 3. Again, the left panel contains all unstable modes, and the right panel, 
one model with logTcff = 4.195 and logg = 4.02. The well-known property 
of this type of diagrams is the zero phase difference and the same amplitude 
ratio for all modes with £ = 1. This is because the light variation of the £ = 1 
mode in SPB models is totaly dominated by the temperature effects. We can 
see an overlapping of domains with different values of £. This can be partly 
removed by considering only one model (the right panel). Another instructive 




9b - % [rad] - 9^ [rad] 

Figure 3: The same as in Fig. 2 but for Slowly Pulsating B-type star models of 5 
Mq. The left panel contains all models in the main sequence evolutionary phase and 
the right panel, the model with logTcfi = 4.195 and logg = 4.02. 

information can be drawn from instability conditions. Fig. 4 shows the amplitude 
ratio and the phase difference as a function of frequency for the model with 
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Figure 4: Amplitude ratios (on the left) and phase differences (on the right) in the BR 
bands as a function of oscillation frequency for the SPB model with logTcff = 4.195 
and logg = 4.02. 

logTcff = 4.195 and logg = 4.02. As one can see, modes with £ = 1 are 
unstable only for the lowest frequencies {ly < 0.9 c/d). The instability is shifted 
to the higher frequencies for higher £-modes, e.g., the £ — 5 modes become 
unstable for v > 1.1 c/d. 

6 Scuti pulsators 

As representatives of S Set pulsator, I took models with M = 2Mq in the 
main sequence phase. The corresponding logTcff range is (3.963, 3.854). All 
calculations were made under assumptions of the mixing length theory and 
the convective flux freezing approximation. The mixing length parameter was 
ctconv = 0.0, i.e., an assumption of inefficient convective transport. The photo- 
metric BR diagram for this type of pulsators is presented in Fig. 5. All unstable 
modes of M = 2Mq models are shown on the left hand side, and modes in 
the (logTcff = 3.909, logg = 4.02) model on the right hand side. As we can 
see, there is some overlap, especially for the £ < 2 modes. Fixing the model 
parameters helps in removing this ambiguity. 

The common opinion about 6 Set pulsation modes is that the main infor- 
mation about the £ values is contained in the phase differences. In fact, this is 
true only for low degree modes with £ < 2, but for higher ^'s we can learn much 
also from the amplitude ratios. In Fig. 6, the amplitude ratio (the left panel) 
and phase difference (the right panel) as a function of oscillation frequency are 
depicted. As we can see, not much can be we achieved for low degree modes 
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Figure 5: The same as in Fig. 2 but for S Scuti star models of 2 Mq. The left panel 
containes all models in the main sequence evolutionary phase, and the right panel, 
the model with logTeff = 3.909 and logg = 4.02. 
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Figure 6: Amplitude ratios (on the left) and phase differences (on the right) in the BR 
bands as a function of oscillation frequency for the <5 Set model with logTcfi = 3.909 
and logg = 4.02. 



also from the Ab/Ab vs. frequency plot because the £ <2 modes are excited 
with very close frequencies. Another interesting property to note is that for 
very high frequency mode {ly > 35 c/d), the phase differences are almost the 
same for all degrees i. 
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Adding radial velocity measurements 

The radial velocity variation averaged over stellar disc is expressed by the well- 
known Dziembowski's (1977) formula 



GM 



where 



and 



h\{fi)fi'^Pi{fi)dfi, 



(7) 
(8a) 

m 



From observations, the radial velocity variations are determined by calculating 
the first moment, A4i, of a well isolated spectral line. 

In this section, I show examples of diagnostic diagrams constructed from 
the radial velocity variation and the light variation in the Johnson R filter. 
Unstable oscillation modes for /3 Cep, SPB and 6 Set star models considered in 
the previous section are presented in Fig 7, 8 and 9, respectively. As we can 
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Figure 7: The location of unstable modes with degrees up to 6 for /? Cep models 
of 12 A/0 on the diagnostic diagrams constructed with amplitudes and phases for 
the R passband and the radial velocity variations. As in Fig. 2, the left panel shows 
all main sequence models, and the right panel, one model with logTcff — 4.400 and 
logs = 3.89. 
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Figure 8: The same as in Fig. 7 but for SPB models of 5 Mq. In the right panel 
unstable modes for the model with logTcff = 4.195 and \ogg = 4.02 are plotted. 
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Figure 9: The same as in Fig. 7 but for 5 Set models of 2 Mq. In the right panel 
unstable modes for the model with logTeff = 3.909 and logg = 4.02 are plotted. 



see, the configuration of £ modes domains differs from that in the photometric 
diagrams. In particular, there are larger phase differences between light and 
radial velocity variation than between two passnands. 

Moreover, with two photometric passbands and the radial velocity data, 
one can apply the method of simultaneous extracting from observations the 
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degree £ and the nonadiabatic parameter /. Determination of the empirical 
value of / allows to avoid, in the process of the ^-identification, the input 
from pulsation models, which still needs many improvements, for example the 
pulsation-convection interaction, opacity tables, effects of diffusion, mixing etc. 
On the other hand, the parameter / constitutes a new asteroseismic probe, 
giving information on subphotospheric layers and is complementary to oscillation 
frequencies determined by stellar interior. Comparing empirical and theoretical 
values of /, one can draw conclusions about the efficiency of convection in 6 
Set stars (Daszynska-Daszkiewicz et al 2003,2005b) or about opacities in j3 Cep 
stars (Daszynska-Daszkiewicz et al. 2005a). 

Uncertainties 

There are many sources of uncertainties in all results presented in this paper. 
Firstly, all calculation were done assuming the zero rotation approximation. Ef- 
fects of rotation can spoil the nice mode separation in the diagnostic diagrams, 
making them dependent on the inclination, rotation rate and the azimuthal 
order m. One such effect is rotational mode coupling. It takes place if the 
frequency distance between modes, with the degrees (. differing by 2 and the 
same to, is of the order of the rotation frequency. In such a situation, the pho- 
tometric amplitude for the coupled mode has to be calculated as superposition 
of all mode amplitudes satisfying the conditions £k = + 2 and ruk = rrij. 
The effect of rotational mode coupling on the photometric diagrams for (3 Cep 
star was studied by Daszynska-Daszkiewicz et al. (2002). Examples for S Set 
stars can be found in Daszynska-Daszkiewicz (2007). Another case when the 
rotation has to be included is when the pulsational frequency, lu, is of the or- 
der of the rotational frequency, so that the perturbation approach fails. It 
happens often in the case of the rapidly rotating SPB stars, in which high or- 
der gravity modes are excited. Such slow modes are no longer described by 
the spherical harmonics, and more complicated formalisms are needed. One 
possibility is the use of the traditional approximation which allows expressing 
the angular dependence of eigenfunctions by the Hough functions (e.g. Lee Si 
Saio 1997, Townsend 2003b). The formula for the light variation due to low 
frequency oscillation was given by Townsend (2003a). Daszynska-Daszkiewicz 
et al. (2007) discussed a prospect for mode identification from diagnostic dia- 
grams and derived the expression for the radial velocity variation. 

Another uncertainty comes from the input physics and atmospheric models. 
The effect of metallicity parameter, Z, and opacities on the diagnostic diagrams 
for P Cep star models was considered by Cugier et al. (1994). For a higher 
value of Z, the separation of modes with different ^'s is much better. Similarly, 
computations with the OP tables, instead of OPAL, lead to a little better £ dis- 
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crimination, especially the radial modes are much more spread. Then, Cugier 
& Daszyriska (2001) checked the effect of the atmospheric metallicity param- 
eter, [m/H], and the microturbulence velocity, ^t, on the diagnostic properties 
of photometric diagrams. For f3 Cep stars these parameters have negligible 
impact. 

The main problem in applying the method of diagnostic diagrams to 6 Set 

modes is that photometric amplitudes and phases exhibit a strong dependence 
on the treatment of convection (Balona &l Evers 1999). To circumvent this 
problem, Daszyriska- Daszkiewicz et al. (2003, 2005b) invented the method of 
simultaneous determination of degree £ and the nonadiabatic parameter / from 
observations. Then we can identify the £-value independently of the pulsation 
models and, by comparing empirical and theoretical / values, constraints on 
convection can be inferred. The result was that the convective transport in S Set 
stars studied by us is rather inefficient. The progress in modelling 6 Set pulsation 
with time-dependent convection treatment was achieved e.g. Grigahcene et al. 
(2005), Dupret et al. (2005a) and Dupret et al. (2005b). 

In the case of 6 Set models, the uncertainties in atmospheric models can 
play much more important role than in the B-type pulsators. In Kurucz stan- 
dard models, the flux derivatives over effective temperature and gravity are 
not smooth at the temperature where convective transport becomes important. 
This is because of using an overshooting approximation that moves the flux 
higher in the atmosphere, above the top of the nominal convection zone. In 
the NOVER-ODFNEW models computed by F. Castelli the problem of non- 
smooth derivatives does not exist. There are also NEMO models (New Model 
Grid of Stellar Atmospheres) which include different treatment of convection. 
The models were computed by the Vienna group with modified versions of the 
Kurucz ATLAS9 code. The grids have smaller steps in Tes and logp than in Ku- 
rucz's computations, and the flux derivatives are perfectly smooth. The effect 
of using various atmospheric models in the calculation of 6 Set observables was 
discussed by Daszyriska-Daszkiewicz et al. (2004) and Daszyriska- Daszkiewicz 
(2007). 

Conclusions 

There is a potential for mode identification from the BRITE photometry. As we 
could see, two-colour information can yield some constraints on the spherical 
harmonic degree, i, of the oscillation modes excited in main sequence pulsators. 

However, space observations should be followed up by simultaneous ground- 
based spectroscopy. The advantages of adding the spectroscopic variations is 
obvious. By combing photometry and the radial velocity data, we will improve 
significantly the discrimination of i and infer better seismic constraints on stel- 
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lar parameters and input physics. Moreover, from the line profile variations the 
identification of the azimuthal order, m, becomes possible. With such unprece- 
dented data, we can hope for a great step In asteroseismology of main sequence 
pulsators with the BRITE-Constellation. 
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